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ABSTRACT 



Filmwise condensation of steam on externally-finned tubes 
is a very complex process. Recent experiments have shown 
that enhancement ratios (ratio of steam-side heat-transfer 
coefficient to that of a smooth tube having the same diameter) 
exceeded the area enhancement produced by the fins. Moreover, 
the enhancement ratios for fully flooded tubes exceeded the 
values predicted by a simple, one-dimensional conduction 
model by a factor of 2 to 4 . A new two-dimensional conduction 
model was developed, which showed that the one-dimensional 
model overpredicted the two-dimensional results for high 
conductivity tube-metals such as copper by as much as 13%. 

The two-dimensional model also showed that variations in fin 
thickness or spacing can result in an overprediction by the 
one-dimensional model of the two-dimensional results by as 
much as twenty- two percent. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Because the size and weight of main propulsion and 
auxiliary systems are inversely proportional to the ease in 
which heat or energy is exchanged, a major thrust of the 
research and development program of the U.S. Navy is toward 
the design and development of equipment that perform at 
higher efficiency with a reduction in size and weight. In 
particular, one area of great interest is the reduction in 
size and weight of steam condensers. 

Current equipment is of the bare or smooth tube type of 
heat exchanger in which the steam is condensed on the outside 
of the tubes (the shell side) . The heat-transfer effective- 
ness is governed by the amount of surface and the overall 
resistance to the flow of heat. This overall resistance, in 
turn, is the sum of the vapor-side, tube wall, and water-side 
thermal resistances. Reduction of any or all of these 
resistances results in an increase in heat-transfer performance 
with an associated reduction in condenser surface area which 
translates into decreased size. 

Both the vapor-side and tube-side thermal resistances can 
be considered as the reciprocal of the so-called hS product 
where h is the coefficient of heat transfer and S is the 
surface area through which the heat flows. Because of the 
curved surface of the exchanger tubing, the vapor-" space" 
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and tube-side surfaces are not equal and an adjustment must 
be made to the tube-side heat-transfer coefficient to allow 
for a design on the basis of the vapor or shell-side surface. 

It is a fact that the vapor-side hs product (reciprocal 
of resistance) can be as high as two or three times the tube- 
side hS product. Thus, it would appear reasonable, in an 
attempt to decrease the overall resistance, to attack and try 
to increase the hS product on the tube side. Unfortunately, 
little can be done in the arena of increasing h significantly 
on either side although the vapor side has greater possibili- 
ties in this regard and, by its very location, not much can 
be done with regard to the increase of the tube-side surface. 

It is for the foregoing reasons that recent efforts have 
focused on the use of externally-finned tubes in steam con- 
densers to increase the flow of heat by increasing the vapor- 
side hS product [Refs. 1,2,3]. 

For the past few decades, it has been accepted that the 
use of externally-finned tubes would be inappropriate mainly 
because of the liquid retention (flooding) that occurs between 
fins (Figure 1.1). This flooding is a result of the relatively 
high surface tension of water (at 50°C water has a surface 
tension of 0.068 N/m, while refrigerants possess values about 
0.015 N/m). One of the unavoidable results of flooding is the 
considerably lower heat transfer through the flooded portion 
than that through the unflooded portion. For the case of an 
externally finned copper tube with steam as the condensate. 
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External Diameter of Fins 




Figure 1 . 1 . Schematic of Condensate Retention on Finned Tubes 
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the copper has a thermal conductivity 600 times greater than 
that of water (385 W/m*K for copper compared with 0.6 W/m-K 
for water) . Due to its overall importance in the heat- 
transfer scheme, the flooding of the lower portion of the tube 
by condensate has been studied theoretically and experimentally 
by several investigators [Refs. 4, 5, 6, 7, 8]. These studies have 
shown that the flooding angle depends mainly on the fin 
spacing and the ratio of surface tension to density of 
condensate as given by the following equation for the case of 
rectangular-section fins [Ref. 5] : 

Ip = cos" 1 [1 ^£_] (1.1) 

r pgDs 

Recent work by Yau et al. [Ref. 9] and Wanniarachchi et al. [Ref. 1] 
has shown that even in the cases with considerable flooding, 
the heat-transfer performance exceeds that expected on the 
basis of area increase due to the fins. 

In 1948, Beatty and Katz [Ref. 10] developed the first model 
to predict the condensation heat-transfer coefficient for 
horizontal, externally-finned tubes. Their model ignores 
the existence of flooding. Although their model agrees to 
within 10 percent for condensing fluids possessing low surface 
tensions, it tends to be rather inaccurate for high-surface- 
tension fluids, such as water. Several other investigators, 

Rudy and Webb [Ref. 2] and Owen et al. [Ref. 8] attempted refinement 
of the Beatty and Katz model but still were unable to predict 
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the heat-transf er coefficient to a reasonable degree of 
accuracy . 

Rudy and Webb used the Beatty and Katz model for the 
unflooded portion, while neglecting heat transfer through the 
flooded portion. With this modification, Rudy and Webb under- 
predicted data for low-surface-tension fluids by up to 
50 percent. Owen et al . attempted to improve the Rudy and 
Webb model by allowing heat transfer through the flooded 
portion. In order to compute heat transfer through the 
flooded portion, they assumed parallel conduction paths 
through the fins and the condensate that is flooded in the 
interfin space as shown in Figure 1.2. With this modification, 
they improved the predictions slightly; underprediction for 
low-surface- tension fluids decreased to about 40 percent. 
However, as shown by Honda and Nozu [Ref. 11], both of these 
models overpredicted the steam data of Yau et al. by a factor 
of up to 3. The fact is that the use of Beatty and Katz 
model, which does not take surface tension forces into consi- 
deration for the unflooded portion, may be responsible for 
this large discrepancy. 

The most recent and accurate attempt at modelling the 
problem was accomplished by Honda and Nozu [Ref. 11] . Their 
resulting solution, a fourth-order differential equation for 
the film thickness, was solved numerically subject to reasonable 
boundary conditions. An agreement to within ±20% was achieved 
by their model in predicting data experimentally obtained 
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Figure 1.2. One-Dimensional Parallel Conduction Kodel 
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for 11 fluids and 22 tubes. However, the disagreement with 
steam condensation data was as high as 40%. 

To understand the condensing heat transfer performance 
more clearly, it is convenient to define the enhancement ratio 
as proposed by Wanniarachchi et al . iRefs. 1,2]. The enhancement 
ratio is defined as the ratio of the steam-side heat-transfer 
coefficient of a finned tube to that of a smooth tube with the 
same diameter as the finned tube root diameter at the same 
heat flux. 

The present investigator has been mainly motivated by the 
heat transfer through a finned tube that is completely flooded 
by condensate. As shown by Winniarachchi et al. [Ref. 2], a finned 
tube with a fin spacing of 0.5 mm, and a height and thickness 
of 1.0 mm and 1.0 mm respectively showed a considerable 
improvement in heat-transfer performance even though this tube 
was completely flooded by condensate. In fact, this tube 
showed an enhancement ratio of about 2.5 under vacuum (85 mmHg) 
and 4.5 at atmospheric pressure. These ratios even exceeded 
the area enhancement (A^/A g ), which was 2.11. If one- 
dimensional conduction is assumed (see Owen [Ref. 8] ) , the expected 
enhancement ratio was computed to be about 0.8. This unex- 
pectedly high enhancement ratio may be, at least in part, 
explained by two-dimensional conduction and convection in 
the interfin space. 
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B. OBJECTIVES 



The main objectives of this thesis are: 

1. Develop a two-dimensional numerical model of a horizontal- 
finned tube fully flooded with steam condensate. The 
program using the model shall be interactive to allow 
different cases of fin thickness, fin spacing, and 
tube-metal thermal conductivity to be studied, 

2. Using a thermal analyzer computer program, determine 
the temperature distribution plots and heat flux for 
each case, 

3. Compute enhancement ratios for each one of the cases 
studied and compare these ratios with those computed 
by using a one-dimensional model, 

4. Compare the one and two-dimensional enhancement ratios 
and comment on the importance of the role of two- 
dimensional heat transfer in the overall heat-transfer 
performance, 

5. Determine the value and usefulness of the two-dimensional 
model as an aid in the prediction of enhancement ratios 
for various cases. 



20 



II. DEVELOPMENT OF THE TKO DIMENSIONAL MODEL 



A. THE PRINCIPLE OF THE FINITE DIFFERENCE 

1. Establishment of the Spatial Increment 

Consider the three-dimensional Laplace steady-state 
equation for the spatial temperature distribution in a 
configuration 






0 



( 2 . 1 ) 



It may be solved analytically, subject to certain boundary 
conditions, and the solution is called an analytical solution. 
The finite-difference solution is arrived at numerically and 
applies to a network of points called nodes, with the 
temperatures at the nodes being estimated. The numerical 
solutions obtained by analytical and finite-difference 
methods can be made to agree with any desired accuracy. 

Consider a grid of nodes. Applications of equation 
(2.1) at a particular node whose temperature is T q allows 
estimation of the second derivatives in terms of the 
temperature T q and the six related temperatures at the 
adjacent nodes of the grid. In Figure 2.1 an arbitrary 
function T = f(x) is shown which considers only the trace 
of T = even though equation (2.1) shows that there 

will also be a functional relationship between T and v, 
namely, T = f^(y) and between t and z, namely, T = f^tz). 
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fe m p e rai ur e. 




Figure 2.1. Approximation of the First Derivative 
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If y and z are held constant, the first derivative 
of equation (2.1) may be approximated by the change in T 
divided by the change in x: 



3T ~ AT + AT 
3x AXj + Ax 2 




- T ) + (T 

o o 

Ax^ + Ax 2 





- T 



1 



+ Ax 2 



( 2 . 2 ) 



An estimate of the second derivative is obtained by 
estimating the first derivative in each of the two intervals 
and dividing the difference of the first derivatives by the 
total interval over which the first derivative estimates 
apply. The total interval will be the distance between the 
centers of the two intervals, Ax^ and Ax 2 . Hence 




AT 2 /Ax 2 - AT 1 /Ax 1 




If A*! = Ax 2 = Ax, 



3 2 T _ 

3x 2 



(T 2 -T q )/ax - (T q -T 1 )/Ax 



Ax 



T. + T_ - 2T 
1 2 O 

Ax 



(2.3) 



A similar development using T = f 2 (y) with temperatures 
T^ and T^ yields 




T 3 + T 4 



- 2T 



ay 



(2.4) 
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and for T = f^(z) with temperatures and 



a 2 T 

9z 2 



T c + T, - 2T 

3 D O 

Az 



(2.5) 



so that equation (2.1) may be approximated by 



T. + T 0 - 2T T-. + T. - 2T T c + T, - 2T 
1 2 o + ^ 4_ o + _J 6_ ° =0 (2.6) 



Ax 



Ay 



Az 



and if Ax = Ay = Az, the result is 



T 1 + T 2 + T 3 + T 4 + T 5 + T g - 6 T q = 0 



Further justification of the result indicated by 
equation (2.3) comes from a consideration of a Taylor Series 
Expansion. Recall that the Taylor Series Expansion is given 
by 



- f (a) .0 ^ f' (a) 



f" (a) 



f (x) = — Qj— (x-a) + — j-j — - (x-a) + — ^ — L (x-a) 



f (a) , . 3 

+ — — -(x-a) + . 



<r n / \ 

, f (a) , , n 

+ -jn— (x-a) 



(2.7) 



Two cases will be examined. First for f (x) = T(x) 



at x = x^ , then 



T(x=x 1 ) = T 



T(x=x ) = T 
o o 
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a = x 



x-a = Ax 



and with these in equation ( 2 . 7 ) 



T, = 



m ^ dT 




i 


d 2 T 


2 


i 


3 

d T 


— T 4- -r 

o dx 


Ax 

X 


+ 2 


d x 2 


Ax 

X 


+ 6 


, 3 
dx 



Ax' 



x 



1 d 4 T 



24 , 4 
dx 



Ax + 



Then for f (x) = T(x) at x = x. 



T(x=x 2 ) = T 2 

T(x=x ) = T 

o o 



a = x 
x-a = -ax 

With these in equation (2.7) 



T 0 = t 
2 o 



dr 

dx 



Ax + IT 



X 



1 d 2 T 



^ dx^ 



Ax 



1 d 3 T 



X 



6 - 3 
dx 



Ax' 



x 



4 

1 d T 



24,4 

dx 



Ax 



x 



When equations (2.8) and (2.9) are added , one 



T 2 + T 1 



. d 2 T,. x 2 , 1 d 4 T /A x 4 , 
2 T q + 2 ^ Ax ^ + 12 — 4 ( Ax ) + 



dx 



dx 



( 2 . 8 ) 



(2.9) 

obtains 

( 2 . 10 ) 
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If Ax is small, then all terms above the second order 



in Eq. (2.9) may be considered negligible. Thus upon 
rearrangement, it is seen that 



d 2 T 

dx 2 



T 1 + T 2 - 2T o 
Ax 2 



( 2 . 11 ) 



which is a result that is identical to eauation (2.3) . Of 
course a similar development will provide a substantiation 
of equations (2.4), (2.5) and (2.6). 

2 . Node to Datum Analysis 

Using the approximations for the second derivative, 
a node analysis can be performed to predict the temperature 
distribution in a mathematical model of a physical configura- 
tion. With a node analysis, the configuration of interest 
is divided into small but finite subvolumes with each sub- 
volume being considered as isothermal. It may be then 
assumed that the geometric center of each subvolume is 
representative of that subvolume and this geometric center 
is called the node-point or node. 

A node analysis consists of an application of the 
First Law of Thermodynamics to generate energy balances for 
each node. This is done for each node by equating 



HEAT IN + HEAT GENERATED = HEAT OUT + HEAT ACCUMULATED (2.12) 
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For the two dimensional case displayed in Figure 2.2, the 
development of the node equations for node point (j,k) is 
shown in Figure 2.3. Assuming heat flow by conduction and 
a nodal dissipation, q, a steady state or First Law energy 
balance requires that the sum of all heat flows (Watts or 
BTU/hr) leaving the node equals zero. Thus with a heat input 
as a negative heat flow outward, 



q x + q 2 + q 3 + q 4 - qi = 0 



(2.13) 



In terms of the temperatures in Figure 2.4 and in terms of the 
Fourier Law 



KA (AT) 
L 



*1 

^2 

^3 




KAx 

Ay 



KAy 

Ax 



6 (T -T, ) 
o 1 

6 (T -TO 
o z 

6 (T -TO 
O 3 



^4 



KAx 

Ay 



6(T o -T 4 ) 



(2.14) 



and when these are inserted into equation (2.13), the result 
is 
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Figure 2.2. Two-Dimensional Case 



28 



o 



(j>+l) 

o 



o 

ti 




0,M) 



Figure 2.3. Node Equation Point (j,k) 
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Figure 2.4. Node Equation 
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Ti 4 <V T i> + T7 6 (t o- t 2> + ^ 6 'W 




(2.15) 



If two conductances are defined by Kl = KAy6/Ax and 
K2 = KAx6/Ay then 



KKT^) + K2(T q -T 2 ) + K1(T q -T 3 ) + K2(T q -T 4 ) = q i 



(2.16) 



and if there is no heat generation (q^=0) , the equation 
becomes 



T q ( 2K1 + 2K2) - KK^) - K2(T 2 ) - Kl (T 3 ) - K2(T 4 > = 0 

(2.17) 

For the case where there is heat flow by both 
conduction and convection as in Figure 2.5, with a heat 
input as a negative flow outward. 



For the conduction terms, q^ , q 2 and q 3 , equation (2.141 is 
still applicable. However, for the convection at the boundary 
as indicated by q c (Figure 2.5) 




0 



(2.18) 



<3, 



hAyS (AT) 



(2.19) 



c 
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Figure 2.5. Node on Convecting Boundary 
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In terms of the temperatures in Figure 2.5 



K (Ax/2) 



* 6 (T -T n ) 

Ay o 1 



*2 m 



q 3 = ^f^6(T 0 -T 3 ) 

q c = MyS|T o -T s ) 

and when these are inserted into equation (2.18), the result 
is 



K(A *^ 2) 6 (T -T ) + 6 (T -T ) + K 6 (T -T.J 

Ay o 1 Ax 0 2 Ay o 3 



+ hAy6 (T q -T s ) = q ± 



( 2 . 20 ) 



Letting Kl = K - ( ^ — , K2 = and K3 = hAy6 , then 



Kl (T -T, ) + K2(T -T ) + Kl (T -T,) + K (T -T ) = q. 

O X U Z U J U o X 



( 2 . 21 ) 



which can be simplified to 



(2K1 + K2 + K3)T q - K1T 1 - K2T 2 - KlT 3 = q ± 



( 2 . 22 ) 
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(2.23) 



and once again if there is no heat generation (q^ =0), the 
equation becomes 

(2K1 + K2 + K3)T - KIT, - K3T- - KIT- - K3T = 0 

o i 2 3 s 

For the case of a node on an adiabatic boundary with no 
heat generation as in Figure 2.6 



+ q 2 + q 3 = o 



(2.24) 



In terms of the temperature in Figure 2.6 and in terms of the 
Fourier Law 



- T'w 

<12 - “W 

<33 “ ^ ‘V3 1 



and when these are inserted into equation (2.24), the result 
is 



TT M V T l ) + ^ 4(t o- t 2> + K(fl Ay 2> S< V T 3> 



0 

(2.25) 



Letting K1 = 



K (Ax/2) 
Ay 



6 



and K2 



KAy 6 
Ax 



(2K1 + K2)T q -KlT^ -K2T 2 - KlT 3 = 0 (2.26) 

It is easy to see that when one proceeds in this 
manner for all nodes, an nxn set of simultaneous algebraic 
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Figure 2.6. Node on Adiabatic Boundary 
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equations will result. These may be represented in matrix 
form as 



[K] [T] = [Q] (2.27), 

where [K] is an nxn matrix of coefficients and IT] and IQ] 
are nxl column vectors representing the node temperatures 
and heat inputs respectively. 

The solution of equation (2.27) is handled by the 
computer program TVSSI [Ref. 12:Appendix A]. 

B. MODEL BUILDER TVIN 

In order to provide input to program TVSSI, it was 
necessary to create a model of the configuration to be 
analyzed. The overall model of the condensing region 
incorporated three major features. Because the system to 
be modelled was an externally finned tube with the external 
space flooded by condensate, the azimuthal portion of the 
model was chosen to be 1/2 tt radians. This is shown in 
Figure 2.7 and because the areas of each node contain a tt 
term, the use of 1/2tt simplified the construction of the node 
equations by permitting cancellation of the it terms. The 
second feature required that the maximum number of nodes to 
be used for the model was 300, due to input restrictions 
(the dimension statement) of the main thermal analyzer. The 
final requirement necessitated that the program be written so 
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Figure 2.7. Model Azimuthal Section 
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as to allow the user to input the variables listed in Table I 
interactively . 

The basic model is shown in Figure 2.8. The upper 
horizontal edge represents the condensate/vapor interface 
while the lower horizontal edge represents the inner wall of 
the tube. Heat-transfer coeffients were specified for the 
inner-tube wall and the vapor-condensate interface. The 
model assumes that the heat is transferred across the retained 
condensate by conduction only. Due to symmetry, the dotted- 
dashed lines at the two vertical edges of the figure represent 
an isolation of a f in-condensate-tube entity into a repeating 
section. These lines, of course, represent adiabatic 
boundaries. Thus, the model assumes that the only heat 
transferred into or out of the model occurs at either the 
condensate/vapor interface or the inner tube wall. 

The first model constructed consisted of 299 nodes with 
13 columns and 23 rows. The outermost rows and columns 
consisted of half node subvolumes. This facilitated the 
placement of the node centers at both the inner tube wall 
and upper condensate/vapor interface in addition to the 
boundary of the outer tube wall and fin perimeter. Nine rows 
were placed in the condensate region above the fin and in the 
region including the fin and adjacent condensate. The 
remaining four rows were placed in the tube wall region. 

In addition the fin region contained 4 columns of nodes 
while the remaining eight columns were placed in the inter-fin 
space . 
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Figure 2.8. Basic Model 



TABLE I 



INPUT VARIABLES FOR MODEL BUILDER TVIN 



T^ - Temperature of fluid inside tube 

T q - Temperature outside condensate layer 

r. - Inner radius of tube 
1 

r - Outer radius of tube 

0 

h. - Heat transfer coef. at inside tube wall 

1 

h Q - Heat transfer coef. at condensate layer 

- Thermal conductivity at tube and fin 

- Thermal conductivity of condensing liquid 

e - Fin height 

t - Fin width 

s - Spacing between fins 

z - Height of condensate above tip of fin 
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The location of these columns and rows determined the 
height and width increments of the various subvolumes, 
depending on their location in the model. For example, if 
the tube wall thickness is four millimeters, and the tube 
wall area contains four rows, the height increment for a 
node in this region would be 4 mm/ 4 rows =lmm. The width 
increments for a node would be computed in identical fashion 
by dividing the size of the particular region by the number 
of columns occupying that region. Figure 2.9 identifies all 
of the different height and width increments utilized by the 
model . 

C. EVALUATION OF RESULTS 

1 . Verification of Accuracy 

Verification of the model consisted of using the 
temperatures generated by TVSSI (Using TVIN as input) and 
the conductance values employed in TVIN to compute the heat 
flow across three boundaries in the model: 

(1) condensate upper surface 

(2) fin perimeter and tube outer wall 

(3) tube inner wall 

The foregoing heat flows were computed in each case 
by using the node temperatures generated by TVSSI with the 
associated internode conductance values used in the input 
generated by TVIN. These heat flows, of course, must be 
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Figure 2.9. Model Increments. (A) Height and (B) Width 



equal and it was observed at the outset that 13x23 model 
yielded certain unacceptable discrepancies in the three values. 

Because the 13x23, 299 node model was found to be 
inaccurate, it was determined that, in order to achieve the 
requisite accuracy, slight adjustment of the number of rows 
and columns was required from run-to-run. 

TVIN was reprogramed to allow flexibility in the 
number of rows and columns utilized in the model as well as 
placement of the nodes in the fin, tube or liquid regions. 

The first change instituted was to fix the number of rows 
contained in the condensate layer above the fin tip at two 
rows. This action was acceptable because the original height 
of this layer was only 0.08 mm compared with a 1.0 mm fin 
height and 3.175 mm tube wall thickness. Fixing the number 
of rows at two permitted the use of additional rows at the 
fin and tube levels. Program TVIN was changed to inter- 
actively accept a specified number of columns (see Appendix 
A) . The number of rows was obtained by making the product 
of rows and columns as close as possible to 300. 

A description of the menu-driven procedure for the 
selection and allocation of rows and columns for a 296-node 
model example now follows: 

If the response to the query 'How many columns are 
desired?' is 8, the program will provide 300/8 = 37.5 or 37 
rows. This yields an 8x37 = 296 node model. 

Suppose the response to the query 'How many rows are 
desired for the tube wall?' is 5. Then the program will 
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allocate 37-5 = 32 rows for the fin and condensate layer 
above the fin. 

However, the number of rows in the condensate layer 
above the fin tip is fixed at 2. Thus the number of rows in 
the fin becomes 32-2 = 30. The result is an 8x37, 296 node model. 
A display of the foregoing model set-up is shown in 
Figure 2.10. 

2 . Model Verifier Ql 

As mentioned in the previous section, once an input 
file had been created by TVIN, and the node temperatures 
were calculated by the thermal analyzer, TVSSI, the solution 
was verified for correctness and accuracy. This check 
consisted merely of a comparison of the heat flow (watts) 
across three surfaces in the model. An error of 3.5% was 
established as the maximum allowable for a solution to be 
considered acceptable. The model verifier Program Ql (see 
Appendix B) , is an offshoot of the model builder TVIN as it 
interactively accepts the same input data required to build 
the model. Once all input parameters are set, the program 
generates and stores each inter-connecting node conductance 
value . 

In the next step, all node temperatures computed by 
the thermal analyzer are read and stored, by location in the 
model. The node numbers of the condensate outer layer, fin 
perimeter and outer tube wall and the inner tube surface are 
determined for each model via an algorithm which requires as 
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Figure 2.10. 8x37, 296 Node Model 
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input, the number of columns and rows being used in the model 
along with their placement. In addition to determining the 
nodes on these surfaces, the algorithm also computes the 
node number of each of the adjoining nodes in the direction 
of the heat flow. Computation of the heat flow by program Ql 
is then simply a matter of taking each surface node conductance 
and multiplying it by the temperature difference between the 
surface node and its adjoining node. The summation of all 
the surface node values results in the total heat flow across 
the particular surface. 

D. THERMAL ANALYZER TVSSI 

The computer program TVSSI used in obtaining the temperature 
distributions was adapted from program TVSS2 listed by Kern 
and Kraus [Ref. 13: Appendix 2a], This adaptation consisted 
of modifying TVSS2 to compute solutions in the SI system and 
accept an input file specifically designed by model builder 
TVIN for each externally finned tube case. This program also 
had to be adapted to the IBM 3033 AP system. 

The program utilizes the Cholesky decomposition scheme, 
and because of the linearization of the radiation terms (not 
used in this study), the program is iterative. A more indepth 
description of the Cholesky Scheme pertaining to the program 
is discussed by O'Hare [Ref. 12:pp. 48-52], 

For all the cases which were studied, the following 
parameters were set for the analyzer: 
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1. The initial temperatures at which all nodes were set 
and at which the computer began the iterative process 
was 100°C. 

2. An eventual accuracy of .05 between the final and next 
to last iterations was used. 

3. A radiation coefficient convergence factor Of 0.66667 
between iterations was used. This is required by the 
program even if radiation is not considered as a mode 
of heat transfer in the model. 

4. The maximum number of iterations that the computer was 
allowed to perform was set at 12 to avoid excessive 
computer running time. 

5. A damping factor of 0.8 was set as an initial damping 
factor based on the number of non-linear terms in all 
of the node equation. 

The ability of TVSSI to correctly solve a numerical model 
was verified by providing as input a data set listed by Kern 
and Kraus [Ref. 13:pp. 326-336] for an extended surface 
problem where the solutions were derived analytically. The 
computer generated solution agreed completely with the 
analytic solution thereby assuring confidence of analyzer 
generated temperatures. 
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III. RESULTS AND DISCUSSION 



A. INTRODUCTION 

During this thesis effort, 15 different cases were studied 
using the programs TVIN and TVSSI. Table II lists the param- 
eters that remained constant in all cases. A summary of all 
cases examined and their results are listed in Table III. 

B. SMOOTH-TUBE EVALUATION 

As a method for comparison of the variation in heat- 
transfer performance, the heat flux for a smooth tube with 
the same diameter as the finned-tube root diameter was estab- 
lished as a baseline value. The ratio of the one- and two- 
dimensional heat fluxes to the heat flux of the finned tube 
with fin width approaching zero were designated and <p ^ 
respectively. The equation used to determine the finned-tube 
of zero fin width heat flow is given bv: 



C. ONE DIMENSIONAL HEAT FLOW 

The one-dimensional model utilized in this study assumed, 
as did Owen et al. [Ref. 8], parallel paths for heat flow across 
the tube wall, fin and retained condensate. Furthermore, a 
simplifying assumption was made that heat flow is purely radial 
and that the heat is transferred across the retained condensate 



Q 



(h 1 ) (D x ) 



D 

o 




(3.1) 
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TABLE II 



CONSTANT PARAMETERS FOR MODEL BUILDING PROGRAM TVIN 



= 30000 W/nrK 

h = 100000 W/m. K 

0 

r . = 6 . 35 mm 

1 

r = 9.525 mm 

o 

e = 1 .0 mm 



= 0.6 W/m- K 

T = 48.0 °C 



T, = 22.0 °C 
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SUMMARY OF CONFIGURATIONS EXAMINED WITH HEAT TRANSFER PERFORMANCE 
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by conduction only. This was accomplished by modifying pro- 
gram TVIN in such a way so that the data file it created con- 
tained only those internode conductance values representing 
heat flow in the radial direction while deleting those conduc- 
tance terms pertaining to axial heat flow. Programs TVSSI and 
Q1 were utilized to compute the node temperatures and heat 
flow for each case. Once again, a maximum error of 3.5% was 
used when heat flows generated by program Ql were verified. 

D. TWO DIMENSIONAL HEAT FLOW 

The heat flow for the two-dimensional case was evaluated by 
the methods discussed in Chapters I and II. 

E. EFFECTS OF FIN SPACING AND FIN WIDTH ON PERFORMANCE 

Heat fluxes were computed for five different cases of fin 
width and spacing. Three cases had a fin width of 1 mm with 
spacings of 0.5, 1.0 and 1.5 mm respectively. The other two 
cases held the spacing constant at 0.5 mm, while the fin width 
was adjusted to 0.5 and 0.25 mm. 

1 . Effects of Fin Spacing 

This section presents results showing the variations 
of heat flux having fin spacing as a parameter. Figures 3.1, 

3.2 and 3.3 present results for tubes having a fin width of 1.0 
mm with fin spacings of 0.5, 1.0 and 1.5 mm. Each figure repre- 
sents the same spacing variation with the tube-metal thermal 
conductivity being varied from 385.0 W/m.K for copper, through 
45.0 W/m K for copper-nickel, down to 16.0 W/m-K for the case of 
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Figure 3.1. Predicted Heat-Transfer Performance of Finned Copper 
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Nickel Tube With 
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Steel Tube With 




stainless steel. 



Similarly, Figures 3.4, 3.5 and 3.6 present 
the ratio of the two-dimensional heat flux to the one-dimensional 
heat flux (Q21) . For each case of tube-metal thermal conduc- 
tivity, the fin spacing of 0.5 mm shows the best heat-transf er 
performance, while the worst performance occurs at a spacing of 
1.5 mm. Figures 3.7, 3.8 and 3.9 present isotherm plots for 
all three cases of thermal conductivity for the finned tube 
have a 0 . 5 mm fin spacing and 1.0 mm fin width. Moreover, the 
ratio Q21 has its smallest value at a find spacing of 0.5 mm 
and increases with increasing fin spacing. Enhancement ratios 
for the various cases had their maximum value at the 0 . 5 mm 
fin spacing and this ratio showed a marked decreaase as the fin 
spacing was increased. 

The results indicate several important facts. The first 
and most important is the significant role the two-dimensional 
heat flow plays in this type of extended surface analysis. The 
increasing Q21 ratio with increasing fin spacing shows that for 
the cases of copper and copper-nickel tube-metal, the axial 
heat flow into the fin upper side surface tends to raise the 

temperature of the fin tip resulting in a choking effect of the 

heat flow into the fin tip surface. This is of enormous conse- 
quence because in each and every case, the majority of the 

heat absorbed by the finned tube is accomplished at the fin 
tip. Because the one-dimensional case allows no axial heat 
flow into the fin side surface, a greater temperature gradient 
occurs at the fin tip resulting in a greater heat flow. As the 



55 



o 



□ 






□ 



□ 




O 




W 



h— ( 

fa 




o 

o 



001 



860 



06*0 _ 88'0 

'4>/ z <P 



080 



840 



56 



Figure 3.4. Comparison of One- and Two-Dimensional Heat-Transfer 
» Performance For Copper Tube With e = 1.0 mm and t = 
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Figure 3.5. Comparison of One- and Two-Dimensional Heat-Transfer 

Performance For Copper-Nickel Tube With e = 1.0 mm and 
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Figure 3.6. Comparison of One- and Two-Dimensional Heat-Transfer Perfor- 
mance For Stainless Steel Tube With e = 1.0 mm and t = 1.0 m 




Figure 3 . 7 . 



Isotherm Plot For Finned Copper Tube 
(Fin and Condensate Regions) 
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Figure 3.8. 



Isotherm Plot For Finned Copper-Nickel Tube 
(Fin and Condensate Regions) 
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Figure 3.9. 



Isotherm Plot For Finned Stainless Steel Tabe 
(Fin and Condensate Regions) 
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fin spacing is increased, the effect of the fin side surface 
is less pronounced since the poor conductivity of the condensate 
inhibits axial heat flow in the region near the fin space 
centerline. This results in a more closer agreement between 
the one-dimensional model and the two-dimensional version. 

This conclusion can be further justified by comparing 
the Q21 ratios for each tube-metal conductivity. The two- 
dimensional effect at the fin tip is decreased as the thermal 
conductivity of the fin material decreases. This does not 
mean, however, that there is no two-dimensional effect for the 
cases with low tube-metal conductivity. As in the case of 
copper, the two-dimensional effect tends to inhibit the heat 
flow into the fin. For the cases of low metal thermal conduc- 
tivity, there is a less pronounced two-dimensional effect at 
the fin tip, but there is a significant two-dimensional heat 
flow into the lower portion of the fin side surface near the 
tube wall. For this reason, the one-dimensional model tends 
to underpredict the two-dimensional model. 

As seen in all the cases for the copper-nickel and 
stainless-steel tube-metal (except the copper-nickel case 
with 0.5 mm fin spacing), all Q21 ratios are greater than 
unity. For the copper-nickel case of 0.5 mm spacing, the reduced 
spacing even with a lower metal conductivity of 45 W/m-K, tends 
to produce a two-dimensional effect at the upper fin tip. 

2 . Effects of Fin Width 

The results of the variation of heat flux with fin 
width as a parameter indicate that fin width plays an important 
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role in the overall heat-transfer mechanism. The two cases 
for each tube-metal thermal conductivity investigated were fin 
widths of 0.5 and 0.25 mm respectively with fin spacing held 
constant at 0.5 mm. For each case of tube-metal thermal con- 
ductivity, the fin width of 0.5 mm displayed one of the highest 
heat fluxes, being second only to the case of the 1.0 mm fin 
width with 0.5 mm fin spacing. The only significant two- 
dimensional fin tip effect can only be observed in the case 
of the copper tube-metal with a fin width of 0.5 mm. 

The same reasoning as in the previous section can be 
applied in this case. Because the fin tip width has been 
reduced, the total amount of heat transferred into the tip 
region is also reduced and the two-dimensional fin tip choking 
effect has less influence. This concurs with the computed Q21 
ratio of 1.04. By reducing the fin width further in each case 
to 0.25 mm, each of the resultant heat fluxes is less than the 
cases with 0.5 mm fin width. For each case of the 0.25 mm fin 
width the Q21 ratio approached 1.3, signifying the comparatively 
large two-dimensional heat flow into the lower fin side surface. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. The two-dimensional numerical model of a horizontal 
finned tube fully flooded with steam condensate developed is 
a useful tool in the study of the relationship between one- 
and two-dimensional heat- transfer effects. 

2. Two-dimensional heat transfer effects play a major 
role in the overall heat transfer mechanism in the case of an 
externally finned tube fully flooded with steam condensate. 

It has been determined that for cases of high tube-metal con- 
ductivity, such as the case of copper, the one-dimensional 
model can overpredict the results of the two-dimensional 
model by as much as 13%, while variation in the fin width or 
spacing can result in the one-dimensional model underpredicting 
the two-dimensional result by as much as twenty-two percent. 

3. Fin width and fin spacing control the magnitude of the 
two-dimensional effect present and the resulting heat- transfer 
performance. An increasing fin width to fin spacing ratio 
will correspond to an increase in the heat-transfer performance. 

4. The one-dimensional model in cases of high tube-metal 
thermal conductivity exhibiting a fin width to spacing ratio 

of greater than unity will overpredict the heat flux predictions 
of the two-dimensional model. 

5. Decreasing the tube-metal thermal conductivity results 
in a degraded heat-transfer performance. 
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6. The two-dimensional fin-tip choking effect is minimal 
in low thermal conductivity tube-metals while the two- 
dimensional heat transfer into the fin lower side surface is 
far more pronounced. 

B. RECOMMENDATIONS 

1. Use the two-dimensional model computed enhancement ratios 
for the various cases to estimate a steam-side heat- transfer 
coefficient and compare with those methods currently being used. 

2. Using the model builder TVIN and thermal analyzer 
TVSSI, conduct a more in-depth study of the effects of variation 
of fin width, height and spacing on the overall heat-transfer 
performance . 

3. Revise the two-dimensional model to include convection 
effects in the retained condensate region. 



65 



noooooooooo 



A£BE5JQJL2£ A 

T VI N COMPUTER GENERATED INPUT ANALYZER PROGRAM 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

PROGRAM TVIN 

PROGRAM TVIN IS DESIGNED TO ACCEPT VARIABLE 
DIMENSIONS AND INPJT PARAMTERS TO CREATE 
AN OUTPUT DATA FILE OF INTER.NODE CONNECTIONS 
AND CONDUCTANCE TERMS. BY ADJUSTING THE 
PRINT STATEMENT, THE OUTPUT WILL EITHER 
REPRESENT A ONE OR TWO-DIMENSIONAL MODEL. 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



C 

c 



DIMENSION CK ! 300 » 4)»T (300) 

INTEGER P ! 300 * 4) ,CO,N,X,ZW, ZW1, ZW2 ,ZW3 ,VI ,V2 ,W,U, 
DOUBLE PRECISION CTHM , CT HL , HI ,HS , RO , R I , A , B , C , D, 

* VAL4.VAL14, VAL6.VAL16, VAL8, VAL1 8 , V AL5 , VAL7 , V AL1 0 , 
*A2»A3,A4,A5»WINC1 , W INC2 , H I NCI , H I NC 2 , HI NC 3 , AL INE , 
*TIN,TS,T,X1 ,X2,VAL2,VAL12,VAL20,A1,RADH, RADL , EL 
CHARACTER *79 TITLE 
CHARACTER *79 FNAME 

DATA IGT,IN, IPR»I WR/6 ,5,4,8/ 

OPEN ( IPR,FILE= 'PRN' ,STATUS= 'NEW' , FOR M= ' FORMATTED' , 

* IOSTAT= ICK) 

IFdCK .NE. 0) WR I TE ( IOT , 20 ) 

'20 FORMAT! • TROUBLE OPENING PRINTER OUTPUT FILE') 
WRITEUOT, 17) 

17 FORMAT!/' ENTER THE TITLE OF THIS STUDY, 79 ', 

* ' COLUMNS MAX. THIS TITLE WILL APPEAR',/, • AT', 

* • THE TOP OF EVERY PRINTED PAGE OF OUTPUT:') 

RE AD (IN, 18) TITLE 

18 FORMAT! A79 ) 

WRITEUOT, 1) 

1 FORMAT!/' ENTER INNER RADIUS OF TUBE:') 

R E AD ( IN , 2 ) RI 

2 F0RMAT(BN,F10. 5) 

WR IT E ! IOT, 3 ) 

3 FORMAT!/' ENTER OUTER RADIUS OF TJ3E : ' ) 

READ ! IN ,2 ) RO 

WRITE! IOT, 4) 

4 FORMAT!/' ENTER DISTANCE FROM OUTER EDGE OF ', 

* 'FIN TO LIQUID SURFACE:') 

READ ( IN ,2 ) A 

WR ITE ( IOT, 5 ) 

5 FORMAT!/' ENTER HEIGHT OF FIN') 

RE AD ! IN , 2 ) B 

WRITE (IOT, 6) 

6 FORMAT!/' ENTER DISTANCE BETWEEN FINS') 

RE AD { IN , 2) C 
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UU MUUU 



WRITE ( IOT, 7) 
FORMAT! / ' ENTER 
RE AD ( IN , 2 ) D 



WIDTH CF FIN') 



WRITE (IOT, 8} 
FORMAT!/' ENTER H 
*' TUBE INNER SURF 
READ (IN, 2 ) HI 

WRITEdOT, 9) 
FORMAT ( / ' ENTER H 
* 'OUTER LIQUID SU 
RE AD ( IN , 2 ) HS 



EAT TRANSFER COEFFICIENT OF 
ACE: ' ) 



EAT TRANSFER COEFFICIENT OF 
RFACE: ' ) 



I 0 
1 1 

12 

13 

14 
1 5 
1 6 

99 



WRITEdOT, 10) 

FORMAT!/’ ENTER THERMAL CONDUCTIVITY OF FINS AND 
*• TUBE:') 

RE AD ( IN , 2 ) CTHM 



WRITEdOT, 11) 
FORMAT!/' ENTER T 
* ' SURROUNDING TJ 
READ ( IN , 2 ) CTHL 



HERMAL CONDUCTIVITY OF FLUID • , 

BE AND F INS: ' ) 



WRITEdOT, 12) 

FORMAT!/' ENTER TEMPERATURE OF FLUID INSIDE TUBE:') 
READ ( IN , 2 ) TIN 



WRITEdOT, 13) 
FORMAT!/' ENTER T 
*• SURROUNDING FLU 
READ! IN, 2J TS 
WRITEdOT, 14) 
FORMAT!/' ENTER T 
♦ • FOR YOUR MODEL: 
READ (IN, 15) N 
FORMAT! BN, 12) 



EMPERATURE EXTERNAL TO ', 
ID:' ) 



HE NUMBER OF COLJMNS DESIRED 

') 



WRITEdOT, 16) 

FORMAT!/' ENTER THE NUMBER OF ROWS DESIRED FOR ', 

*' THE TUBE:') 

READ! IN, 15)U 
WRITEdOT, 99) 

FORMAT!/' ENTER THE NUMBER OF CDLUMNS OESIRED FOR ', 
* ' THE FIN: •) 

RE AD (IN, 15) ZW 



T ( N) ) 

D CK ARRAYS 



CO = INT(300.0/FLOA 
INITIALIZE THE P AN 

DO 25 I=1,C0*N 
DO 25 J = 1 » 4 
P( I, J) = 0 
CK (I , J)=0.0 
CONT INUE 



VALUES COMMON TO ALL RUNS 

L2 1=C0*N 
L 22= 2 
L 2 3= 0 
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L24=0 
L2 5=0 
L26=0 
L27=0 
L2 8= 0 
C 

L3 1=300 
L 3 2= 5 0 
L33=6 
L34=2 
L3 5=4 
L36= 6 
C 

FL41=.01 
FL42=. 66667 
IL43 = 12 
FL44 = .8 
F L 45= 10 0. 

C 

LI NT 1=0 
LINT 2=0 
L I NT 3=0 
C 

C COMPUTE INTERNODE CONNECTIONS 

C 

DO 55 I = 1.C0*N 
X = N* C CO-1 ) + 1 

I F (A MOD (FLOAT! 1-1 ) , FLOAT (N) I.EQ.O) GO TO 30 
IFCAMODCFLOATC I ) , FLOAT (N) ) . EQ .0 ) GO TO 40 
IFCI.LE.N) GO TO 42 
IFCI.GE.X) GO TO 48 
C 
C 

PCI,1)=CI-N)*10+1 
P( I, 2)=CI-1)*10+1 
p ( i, 3)=( i+i mo+i 
P( I » 4 )= ( I +N) *1 0+ L 

GO TO 55 
C 
C 

30 IF(I.LE.N) GO TO 45 
IF(I.GE.X) GO TO 49 
P< It 2 )= ( 1 + 1 ) *10+1 
35 P( I, 1)=( I-N)*10+l 

P( It 3) = (I+N)*10+1 
GO TO 55 
C 
C 

40 IF(I.GE.X) GO TO 51 

IF(I.LE.N) GO TO 47 
P< I » 2 )= ( I- 1 ) *1 0+1 
GO TO 35 
C 
C 

42 PC If 1)=( I-l)*10+l 

P ( If 2)=< I + l)*10+l 
PC If 3)=( I + N)*10+l 
PC 1,41=3014 
GO TO 55 
C 
C 

45 P C I, 1 )= ( 1*20) + 1 

46 PC I, 2 ) = C I + N)*10+l 
PC I, 3) =3014 

GO TO 55 
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c 

c 

47 
C 

c 

48 



C 

C 

49 

50 



C 

C 

51 

C 

C 

55 

C 

C 

C 

C 

c 



c 

c 

c 



c 

100 



c 

105 



C 

110 



c 

115 



p ( i. n=< i-i )*i o+i 

GO TO 46 



P ( It 1 )=( i-Nmo+i 

P ( It 2 )=( 1-1 )*10 + 1 
P ( 1 1 3 )= ( I +1 ) *1 0+1 
P( It 4)=3024 
GO TO 55 



Pi I,2)=(I-1 i*10+l 
P ( It 1 )=( I-Ni*10+1 
P ( It 3 ) = 3024 
GO TO 55 



P( 1,2 )=( 1-1 1*10+1 

GO TO 50 



CONT INUE 

CALCULATE HEIGHT AND WIDTH INCREMENTS 

ZW=INT(FL0AT(N)/3.0) 

Z W 1=N— ZW-1 

Z W2= I NT ((FLOAT (CO ) -10 . 0) /2. 0 ) 

Z W2= 2 

Z W 3= CO- ( U+l )-ZW2 
WINC1=D / FLOAT ( ZWi 
W I NC 2=C/ FLOAT ( ZWI ) 

H I NC 1=( RO— R I )/FLOAT(U) 

H I NC 2 =B / F L 0 AT ( ZW3 ) 

H I NC3=A/ F LO AT { ZW2 ) 

EL = ( WINC1+WINC2J / 2.0 

COMPUTE THE RADIUS OF EACH ROW OF NODES 

DO 1000 1=1, CO 
IFd.GT.li GO TO 100 
R A DH =R I +0 . 5 *H I NCI 
RADL =RI 
GO TO 150 

IF (I .GE.(U + 1) ) GO TO 105 
RADH=RI+( FLOAT ( I ) -0 . 5 ) *H I NC 1 
RADL=RI + ( FLOAT ( I ) - 1 . 5 1 *H I NC 1 
GO TO 150 

IF(I.NE.(U+li ) GO TO 110 
R ADH=RO+H INC2 / 2 . 0 
R ADL=RI + (FLOAT (I )-l«5J*HINCl 
GO TO 150 

IF(I.GE.(U+ZW3iJ GO TO 115 
RaDH=RO+( FLOAT (I i— (FLOAT (U)+0.5) )*HINC2 
R AOL=RO+ (FLOAT (I )- (FLOAT (U) +1.5) )*HINC2 
GO TO 150 

IF (I .NE. (U + 1 + ZW3) ) GO TO 120 
RA0H = R0 + B + HINC3/2 .0 
R ADL = RO+8— H IN C2/2 .0 
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GO TO 150 
C 

120 IF(I.EQ.CO) GO TO 125 

RADH=RQ*B+( FLOAT ( I )-( FLO AT ( U) +0. 5 J -FLOAT (Zrf3 ) )*HINC3 
RADL=RO+B+( FLOAT ( I )-( FLO AT ( U ) +1 . 5 ) -F LO AT { ZW3 ) ) *HI NC 3 
GO TO 150 
C 

125 RADH=RG+B+A 

RADL=RO+B+A-HINC3/2.0 

C 

C 

c 

150 AL INE=RADH**2-RADL**2 

c 

C COMPUTE INTER-NODE CONDUCTANCE \/ALJES 

C 

c 

185 I F ( I . GT . 1 ) GO TO 187 

H T C= H I 
RD=R I 
GO TO 189 

187 IF (I .NE.CO) GO TO 190 
HTC= HS 
RD=RO+A+E 

189 Al=HTC*WINCl*RD/2 .0 
A2=A 1*2-0 
A3=A2*EL/lrtINC 1 
A4=HTC*WINC2*RD 
A5=A4/2 .0 

C 

C 

190 IF(I -GE. (U+2I ) GO TO 200 
C 

VA L4 = RADL* W INC 1*C T HM/ H INC 1 
GO TO 300 

200 IF<I.GE.((U+2)+Zrf3)) GO TO 250 

VAL4=RADL*W INC 1*C T HM/ H I N C2 
GO TO 300 

250 VAL4=CTHL*RADL*WI NC1/HINC3 

300 VAL2=VAL4/2.0 

C 

C 

C 

C 

I F ( I -GE • ( U+l J ) GO TO 310 
VAL14=CTHM*RADH*W I NCI/HI NCI 
GO TO 400 

310 I F ( I . LE • ( U+ZW3 ) ) GO TO 350 
VAL14=CTHL*RADH*H INC1/HINC3 
GO TO 400 

350 VAL14=CTHM*RADH*W INC1/HI NC2 
400 VAL12=VAL14/2.0 

C 
C 
C 

i F ( I .GE • ( U+2 ) ) GO TO 413 

V AL6=CTHM*RADL*EL /HINC1 
GO TO 50C 

410 IF(I.GE.(( U+2 ) +ZW 3 ) ) GO TO 450 

V AL6= ( ( CTHM*W I NC 1 +CTHL*W INC2)*R4DL)/ (2.0*HI^C21 
GO TO 500 

450 VAL6=CTHL*RADL*EL/HINC3 
C 

c 

c 
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500 



61 0 



620 

C 

C 

C 

65 0 



670 



680 

700 

C 

C 

C 



730 



74 0 

750 

C 

C 

C 



76 0 



770 

C 

C 

c 

800 



31 0 



32 0 
C 
C 
C 

845 



850 

C 

C 

86 0 



C 



IF(I.GE.(U+1J i GO TO 613 
VAL16=CTHM*EL*RADH/HINC1 
GO TO 650 

IF ( I .LE . ( U+Zw'3 ) ) GO TO 620 
VAL16=CTHL*EL*RADH/HINC3 
GO TO 650 

V A LI 6 = ( <CTHM*WINC l +CTHL* W IN C2 ) *1 ADH )/ (2.0*HINC2 ) 



IF (I .GE.(U+2i i GO TO 670 
VAL8=RADL*W INC2*C THM/HINCl 
GO TO 700 

IF(I.GE.((U+2)+ZW3)) GO TO 680 
VAL8=CTHL*RADL*WINC2/HINC2 
GO TO 700 

V AL8 =CT HL*RADL *WI NC2/HINC3 
VAL10=VAL8/2. 0 



I F ( I .GE. (U+1J ) GO TO 730 
VAL18=CTHM*RADH*W INC2/HINC1 
GO TO 750 

I F ( I .LE . ( U +ZW3 ) i GO TO 740 
VAL18=CTHL*RADH*W INC2/HINC3 
GO TO 750 

VAL18=CTHL*RADH*W INC2/HI NC2 
VAL2 0=VAL18/2 .0 



I F < I .N£.1(U+1)+ZW3)1 GO TO 760 
VAL5=(CTHL*(RADH**2-(R0+B )**2)+CTHM* ( ( RO+B) 
***2-RADL**2i ) / (2. 0*WINC1 ) 

GO TO 800 

IF(I.GT.((U+1)+ZW3)) GO TO 770 
VAL5=CTHM*ALINE/l 2.0*WINCl) 

GO TO 800 

VAL5=CTHL*ALlNE/( 2 .0*W I M C 1 ) 



I F ( I ,GE-(U+1) ) GO TO 810 
VAL7=CTHM*ALINE/( 2.0*WINC2) 

GO TO 845 

IF(I.NE.(U+1) ) GO TO 823 

VAL7=(CTHL*(RADH* *2— R 0**2 ) + CTHM* ( R0**2-R ADL * *2) ) 
*/ (2.0*WINC2 ) 

GO TO 845 

VAL7=CTHL*ALINE/< 2.0*HINC2) 

COUPLE CONDUCTANCE VALUES TO CONNECTION TERMS 

IF(I.EQ.l) GO TO 905 
IF(I.EQ.CO) GO TO 912 
Q=N* (I — 1 ) 

DO 900 J=Q+1 »Q+N 



IF (J .EU. (Q + l ) ) GO 
I F ( J .LT.( Q+Zri+1) ) 
I F (J .EQ. (Q+Za + 1) ) 
IFIJ.NE. (Q+N) ) GO 



TO 875 
GO TO 880 
GO TO 885 
TO 890 
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c 

C K ( J * 1) = VAL 13 
CK ( J y 21 =V AL7 
C K ( J y 31 =VAL20 
GO TO 900 
C 

c 

880 CK ( J y 11 =VAL4 
CK ( J » 2) =VAL 5 
CK< J,3)=VAL5 
CK< J,41=VAL14 
GO TO 900 
C 
C 

885 CK(J,11=VAL6 
CK ( J » 2 ) =VAL5 
CK ( J y 31 =VAL7 
C K ( J » 4) = VAL 16 
GO TO 900 
C 
C 

390 C K ( J y 1 1 = VAL8 

C K ( J » 2) =VAL7 
CK ( J f 31=VAL7 
CK ( J y 4) =VAL 18 
GO TO 900 
C 
C 

875 CK ( J» 1} =VAL2 

C K { J » 21 =VAL 5 
C K ( J y 3) =V AL 12 
C 

900 CONTINUE 

901 GO TO 1000 
C 

C 

905 DO 910 J=1,N 
IF(J.GT.l) GO TO 906 
C K ( J # 11 =VAL5 

CK ( J y 2 ) =V AL 12 
C K ( J » 3) =A 1 
GO TO 910 
C 
C 

906 IF(J.GE.ZW+1J GO TO 907 
CK (J, 1)=VAL5 

CK ( J y 2) = VAL5 
CK ( J» 31 =VAL 14 
CK(J,41=A2 
GO TO 910 
C 
C 

907 IF(J.NE.(ZW+1) 1 GO TO 908 
CK ( J y 1 1 =VAL5 

CK (J y 21 =VAL7 
CK ( J y 31 =VAL 16 
CK ( Jy 4) =A3 
GO TO 910 
C 

c 

908 IF(J.EQ.N) GO TO 909 
CK(J,11=VAL7 

CK( J,2)=VAL7 
CKlJy31=VAL18 
C K ( J » 4) =A4 
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GO TO 910 
C 
C 

909 CK(J» 1 1 = VAL 7 
C K ( J , 2) = V AL 20 
CK < J , 3) =A 5 

C 

910 CONTINUE 
C 

C 

91 1 GO TO 1000 
C 

C 

912 Q=(C0-1J*N 

DO 917 J=Q+1 « Q+N 
I F C J . GT . ( Q + 1 ) ) GO TO 913 
CK(Jt 1)=VAL2 
CKl J,2)=VAL5 
C K ( J , 3) =A 1 
GO TO 917 
C 

c 

913 IF(J .GE.(ZW+0 + iI) GO TO 914 
C K ( J * 1 J = VAL4 

CK ( J ♦ 2J =VAL5 
C K ( J t 3) =VAL5 
C K ( J » 41 =A2 
GO TO 917 
C 
C 

914 I F C J .NE. (ZW+Q+1) ) GO TO 915 
CK(J, II =V AL6 

C K ( J ♦ 2 J =V AL5 
CM J,3)=VAL7 
C K ( J , 4 ) =A3 
GO TO 917 
C 
C 

915 IFCJ.EQ.IQ + N) ) GO TO 916 
CK (J, 1) =VAL8 

CK ( J t 2) = VAL7 
C K ( J » 31 =VAL7 
CK ( J « 4J =A4 
GO TO 917 
C 
C 

916 CM J , 1) =VAL10 
CK ( J t 2 ) =VAL 7 
CM J,3) =A5 

917 CONTINUE 
10 00 CONTINUE 
C 

C 

C END OF CONDUCTANCE CALCULATIONS 
WRITE (IOT, 922 ) 

922 FORMAT!/, • ENTER NAME OF INPUT FILE':) 

READ! IN, 923 JFNAME 

923 FORM AT( A12 ) 

C 

c 

C CREATE DATA FILE FOR TVSSI 
C 

OPEN (IWR,FILE=FNAME,$TATUS= 'NEW' , FOR M= ' F ORMA T T 
*IOSTAT=ICK) 

C 
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c 

WR ITE (IOT » 925 ) 

925 FORMAT! / » WRITING FILE TVSSI'J 
C 

WRITEdOT *926) CO*N , N» CO, ZW2 

926 FORMAT!/' THIS ',13,' NODE MODEL HAS ',12,' ROWS • 
*» AND ',12,' COLUMNS. THERE ARE ',12,' ROWS') 

WRITE(IOT»927)ZW3 

927 FORMAT!* IN THE CONDESATE LAYER AND ',12,* ROWS ', 
*• IN THE FIN. THE FIN CONTAINS') 

WRITEdOT, 928) ZW 
- 92 8 FORMAT! 12,' COLUMNS. •) 

C 

C WRITE STATEMENTS 
C 

WRITEII WR, 940 ) TITLE 

WRITE (I WR, 944) L21 , L22 , L23 ,L 24,L25, L26 , L27, L28 
WRITE (I WR,950) LI N T 1 , L I NT 2 ,L I NT3 
WR ITE ( I WR, 946)L31 ,L32,L33 ,L34,L35, L36 
WRITE (I WR,947) FL41,FL42, I L43, FL44, FL45 
WR ITE (I WR , 943 ) TI , TS 
C 

DO 1100 1 = 1 »C 0*N 
IF (P!I,4).EQ.O) GO TO 1011 

IF!CK(I,1).GE. 100.0. OR. CKd ,2 ) .GE. 100. 0) GO TO 1005 
IF (CK(1 ,3) .GE. 100 .O.OR.CKII ,4) .GE. 100.0JG0 TO 1005 
WRITE !I WR,948)P! I , 1) ,P( I ,2) ,P! 1 , 3) ,P ! I ,4) 

WRITE! IWR,960) CK( I , 1 ) , CK ! 1 , 2) , CK ! I ,3), CKd, 4) 

GO TO 1100 

10 05 WRITE(IWR,948)P( I ,1),P(I ,2) ,P(I,3) ,P(I,4) 

W R ITE ! I WR , 942 ) CK! 1 , 1 ) , CK ( 1 , 2) , CK ! I , 3 ) , CK ( 1 , 4 ) 

GO TO 1100 

1011 I F!CK! I ,1) .GE .100 .O.OR.CK! I ,2) .GE.100.0) GO TO 1015 
IF!CK!I,3) .GE. 103 .0) GO TO 1015 
WRITEIIWR, 945 )P(I»1),P(I ,2) ,P!I,3) 

WRITE !IWR, 949) CK! I,1),CK(I,2),CK(I,3) 

GO TO 1100 

10 15 WRITE(IWR,945) P( I ,1),P!I,2),P!I,3) 

WRITEIIWR, 983) CK! I ,1) ,CK (I, 2) ,CK(I ,3) 

C 

1100 CONTINUE 
C 

C FORMAT STATEMENTS 
C 

940 FORMATdX, A79) 

941 FORMAT! 6F8. 5) 

942 F0PMAT14F8.4) 

943 FORMAT! 2F8.0) 

C 

944 F0RMAT19I4) 

945 F0RMAT(9I8) 

946 FORMAT ! 9 14) 

947 F0RMAT!F8.3,F8.5, I8,2F8. 1) 

948 F0RMAT!9I8) 

949 F0RMAT!F8.6,F8.5, F8.6) 

950 F0RMATI9I4) 

960 F0RMAT!F8.6,2F8.5,F8.6) 

980 FO RM AT ! 3F8 . 4 ) 

C 

C 

CLOSE(IWR,IOSTAT=ICK) 

IFdCK .NE. 3) WRITE! IOT ,921)EEX 
921 FORMAT! 'TROUBLE CLOSING PRINTER OUTPUT FILE*) 

C 

C 
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STOP 

ENO 
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AEEEMCLiX. fi 

Q1 COMPUTER GENERATED MODEL VERIFIER PROGRAM 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

PROGRAM Q 1 

PROGRAM Q 1 IS DESIGNED TO ACCEPT VARIABLE 
DIMENSIONS, INPUT A DATA FILE OF NODE 
TEMPERATURES GENERATED BY TVSSI AND CREATE 
AN OUPUT FILE OF HEAT FLOW VALUES ACROSS 
I) VAPOR CONDENSATE INTERFACE, 2) FIN AND 
TUBE OUTER PERIMETER, AND TUBE INNER WALL 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



c 

c 



20 



DIMENSION CK( 300, 4) »T (300) 

INTEGER P ( 300 » 4) , CO,N , X, ZW , ZWi , ZW2 , Z W3 , VI , V2 ,W,U, 

DOUBLE precision cthm.cthl, hi,hs ,ro,ri ,a,b,c,d, 

*VAL4, VAL14, VAL6.VAL16, VAL8, VAL18 , V AL 5 , V AL7 , V AL1 0 , 
*A2,A3,A4» A5.WINCI , W INC2 , H I NCI ,H I NC 2 , HI NC 3 , AL I NE , 
*TIN»TS»T»X1«X2,VAL2»VAL12»VAL20»A1,RADH» RADL , EL 
CHARACTER *79 TITLE 
CHARACTER *79 FNAME 

DATA IOT, IN, IPR, I WR/6 ,5, 4,8/ 

OPEN!IPR,FILE=*PRN» ,STATUS= 'NEW' , FORM= • FORMATTED* , 
* I OST AT= ICK ) 

IF ( I CK .NE. 0) WRITE! IOT, 20 ) 

FORMAT! • TROUBLE OPENING PRINTER OUTPUT FILE*) 



WRITEdOT, 17) 

17 FORMAT!/* ENTER THE TITLE OF THIS STUDY, 79 *, 

* * COLUMNS MAX. THIS TITLE WILL APPEAR*, /,* AT', 

* • THE TOP OF EVERY PRINTED PAGE OF OUTPUT:*) 

RE AD ! IN, 18 ) TITLE 

18 FORMAT! A79) 

WRITEUOT, 1) 

1 FORMAT!/' ENTER INNER RADIUS OF TUBE:*) 

R E AD ! IN , 2 ) RI 

2 FORMAT! BN, F10. 5) 

WRITEIIOT, 3) 

3 FORMAT!/* ENTER OUTER RADIUS OF TUBE:*) 

READ (IN, 2) RO 

WRITEUOT, 4) 

4 FORMAT!/' ENTER DISTANCE FROM OUTER EDGE OF ', 

* 'FIN TO LIQUID SURFACE:') 

READ ( IN ,2 ) A 

WRITEUOT, 5) 

5 FORMAT!/* ENTER HEIGHT OF FIN') 

READ I IN, 2 ) 8 

WRITEUOT, 6) 

6 FORMAT!/* ENTER DISTANCE BETWEEN FINS') 

R E AD ( IN , 2 ) C 
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ooon no 



c 

7 
C 

8 
C 

9 

C 

1 0 

c 

1 1 
c 

12 

C 

1 3 

14 
1 5 
1 6 

99 



C 



5 



WRITE! IOT, 7) 

FORMAT*/' ENTER WIDTH OF FIN') 

READ ( IN * 2 ) D 

WRITEdOT, 8) 

FORMAT* / ' ENTER HEAT TRANSFER COEFFICIENT 3F 

* ' TUBE INNER SURFACE: • ) 

RE AD ( IN * 2 ) HI 

WRITE(I0T»9) 

FORM ATI / ' ENTER HEAT TRANSFER COEFFICIENT OF 

* 'OUTER LIQUID SURFACE:') 

RE AD ( IN * 2 ) HS 

WRITEdOT, 10) 

FORMAT!/' ENTER THERMAL CONDUCTIVITY OF FINS 
*• TUBE:') 

R E AD ! IN » 2 ) CTHM 



AND 



WRITEdOT, 11) 

FORMAT!/' ENTER THERMAL CONDUCTIVITY OF FLUID • , 

* • SURROUNDING TUBE AND FINS:') 

READ ( IN * 2 ) CTHL 

WRITEdOT, 12) 

FORMAT!/' ENTER TEMPERATURE OF FLUID INSIDE TUBE:') 
R EAO ( IN, 2 ) TIN 

WRITEdOT, 13) 

FORMAT! /• ENTER TEMPERATURE EXTERNAL TO ', 

*• SURROUNDING FLUID:') 

R E AD ! IN , 2 ) TS 
WR IT E ( I OT , 14) 

FORMAT!/' ENTER THE NUMBER OF COLUMNS DESIRED ', 

* • FOR YOUR MODEL: • ) 

RE AD (IN, 15) N 
FORMAT! BN, 12) 

WRITEdOT, 16) 

FORMAT!/' ENTER THE NUMBER OF ROWS DESIRED FOR ', 

*• THE TUBE:* ) 

READ! IN, 15)U 
WRITEdOT, 99) 

FORMAT!/’ ENTER THE NUMBER OF COLUMNS DESIRED FOR • 
*• THE FIN: •) 

READ ( IN , 1 5 ) ZW 

CO=I NT! 300.0/FL0AT !N) ) 

INITIALIZE THE P AND CK ARRAYS 

DO 25 I=1,CG*N 
DO 25 J = 1 , 4 
P< I, J) = 0 
C K ( I , J)=0.0 
CONTINUE 

VALUES COMMON TO ALL RUNS 

L21=C0*N 
L 2 2= 2 
L23=0 
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L24= 0 
L2 5=0 
L 26= 0 
L27=0 
L2 8= 0 
C 

L3 1=300 
L3 2= 50 
L33=6 
L34=2 
L3 5=4 
L36= 6 
C 

FL41=.01 
FL42=. 66667 
I L43 = 12 
F L44= .8 
F L45= 100. 

C 

LI NT 1=0 
L I NT 2=0 
LI NT 3=0 
C 

C COMPUTE INTERNODE CONNECTIONS 

C 

DO 55 1 = 1, CON 
X=N*(C0-1)+1 

IF(AMOD(FLOAT( 1-1 ) , FLOAT (N) J.EQ.O) GO TO 30 
IF(AMOD(FLOAT ( I ) , FLOAT ( N ) ) . EQ-0 ) GO TO 40 
IF(I.LE.N) GO TO 42 
IF(I.GE.X) GO TO 48 
C 
C 

P ( It 1 1= ( I-N) *10+1 
P ( 1,21=11-1 )*10+1 
P ( I,3)=( 1+1 )*10+1 
P< I, 4)=< I+N)*10+l 

GO TO 55 
C 

c 

30 IF(I.LE.N) GO TO 45 

IF(I.GE.X) GO TO 49 
P(I,2)=(I+1)*10+1 
35 P< I, 1)=(I-N)*10+1 

P ( I, 3 >=( I + N)*10+l 
GO TO 55 
C 
C 

40 IF(I.GE.X) GO TO 51 

IF(I.LE.N) GO TO 47 
PU»2)=(I-1)*10+1 
GO TO 35 
C 
C 

42 P(I, 1)=( I-l)*10+l 

P( I,2i=(I + l)*10+l 
P(I,3)=(I + N)*10+1 
P( I,4)=3014 
GO TO 55 
C 
C 

45 P ( It 1 )= { 1*20) +1 

46 P( I,2)=(I+N)*10+1 
P ( It 3) = 3014 

GO TO 55 
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c 

c 

47 P { I» 1 i= ( 1-1 )*10+1 
GO TO 46 

C 

C 

48 P(I,1)=(I-N)*10+1 
P ( I» 2 ) = ( 1-1 )*10+1 
P ( I * 3 1 = ( I + 1 )*10+1 
P< 1,41=3024 

GO TO 55 
C 
C 

49 P(I,2)=(I-1)*10+1 

50 P ( I, 1 ) = ( I-N)*10+l 
P( I, 3 ) = 3024 

GO TO 55 
C 
C 

51 P( I,2)=(I-1)*10+1 
GO TO 50 

C 

C 

55 CONTINUE 

C 

C CALCULATE HEIGHT AND WIDTH INCREMENTS 

C 

C ZW = I NTIFLOAT (N )/3 . 0) 

Z W 1=N-ZW-1 

C ZW2=INT ((FLOAT (C0)-10.0)/2.0) 

Z W2=2 

ZW3=C0-(U+1)-ZW2 
W INC1=D/FL0AT( ZW) 

W I NC2=C/FL0AT ( ZWI ) 

H INC 1=( RO-RI )/ FLOAT! U) 

H I NC2=B/ FLOAT ( ZW3 ) 

H I NC3=A/F LO AT ( ZW2 I 
EL=( WINC1 + WINC21/2.0 
C 

C COMPUTE THE RADIUS OF EACH ROW OF NODES 

C 

00 1000 1=1, CO 
IF(I.GT.l) GO TO 100 
R A CH=R I +0 . 5 *H I NC 1 
RADL =RI 
GO TO 150 
C 

100 I F ( I .GE.(U+i)J GO TO 105 

RADH=RI +( FLOAT ( I) -0 .5 ) *H INC 1 
R ADL=RI + (FLOAT (I)— 1.5)*HINC1 
GO TO 150 
C 

105 IFU.NE.(U+1) ) GO TO 110 
RADH=R0+HINC2/2.0 
RADL=RI+( FLOAT (I ) -1.5 )*HINC 1 
GO TO 150 
C 

110 I F ( I .GE.U1+ZW3) ) GO TO 115 

RADH=RO + ( FLOAT ( I ) -(FLOAT <U)+0.5) )*HINC2 
RADL = RO+( FLOAT ( I ) - ( FLOAT ( U ) +1.5) J*HI NC2 
GO TO 150 
C 

115 IF(I .NE. (U + 1+ZW3) ) GO TO 120 
RADH=R0+B+HINC3/2 .0 
RADL = R0+B-HINC2/2 .0 
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GO TO 150 
C 

120 IF(I.EQ.CO) GO TO 125 

RADH=RO + B+( FLOAT* I )-( FLO AT < U ) +0 . 5 ) -F LOAT ( ZW3 ) )*HINC3 
RADL = RO+B + (FLOAT( I )-( FLOAT* U ) + 1 . 5 ) -F LOAT * ZW3 ) )*HINC3 
GO TO 150 
C 

125 RADH = RO + B+A 

RADL=R0+B+A-HINC3/2.0 

C 

C 

C 

150 AL INE=RADH**2-RADL**2 
C 

C COMPUTE INTER-NODE CONDUCTANCE VALUES 

C 

C 

185 IFII.GT.1) GO TO 187 
H T C= H I 
RD=R I 
GO TO 189 

187 IF(I.NE.CO) GO TO 190 
HTC=HS 
R D=RO+A+B 

189 Al=HTC*WINCl*RD/2.0 
A2=A 1*2 .C 
A3=A2*EL/HINC1 
A4=HTC*WINC2*RD 
A5=A4/2 .0 

C 

C 

190 I F < I .GE . * U+2 ) ) GO TO 200 
C 

V AL4=RADL*W INC 1*CTHM/HINC1 
GO TO 3 OG 

200 IF(I.GE.((U+2)+ZW3i) GO TO 250 

VAL4=RADL*W INC 1*C T HM/ HINC2 
GO TO 300 

250 VAL4=CTHL*RADL*WINC1/HINC3 
300 VAL2=VAL4/2.0 

C 
C 
C 

c 

IF (I .GE.*U+1) * GO TO 310 
V AL14=C THM*RADH*W INC1/HI NCI 
GO TO 400 

310 I F ( I . LE . ( U+ZW3 ) ) GO TO 350 
VAL1 4=C THL*RADH*W INC1/HI NC3 
GO TO 400 

350 V A LI 4=C THM*RADH*W INC1/HINC2 

400 VAL12=VAL14/2 .0 

C 
C 
C 

I F (I .GE . * U+2 ) ) GO TO 410 
V AL6 =CTHM*R ADL*EL / H INC1 
GO TO 500 

410 IFU.GE.* (U+2J+ZW3)) GO TO 450 

V AL6 = ( * CTHM*W I NC 1 +CTHL*W I NC2 ) *RADL1 / *2. 0*HINC2J 
GO TO 500 

450 V AL6=CTHL*RADL *EL / H INC3 

C 
C 
C 
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500 

61 0 

620 

C 

c 

c 

650 

670 

680 

700 

C 

C 

C 

730 

74 0 

75 0 

C 

c 

c 

76 0 

770 

C 

C 

C 

800 

81 0 

820 

C 

C 

C 

845 

850 

C 

C 

86 0 

C 



IFII .GE. ( U + l) ) GO TO 610 
VAL16=CTHM*EL*RADH/HINC1 
GO TO 650 

I F ( I • LE • ( U+ZW3 ) ) GO TO 620 
VAL16=CTHL*EL*RADH/HINC3 
GO TO 650 

VA LI 6 = { (CTHM*WINC1+CTHL*WINC2)*RADH )/ (2.0*HINC2) 



I F ( I .GE. (U+2) ) GO TO 670 
VAL8=RADL*W INC2*C T HM/ H INC 1 
GO TO 700 

IFd.GE.I (U+2J+ZW3)) GO TO 680 
VAL8=CTHL*RADL*WINC2/HINC2 
GO TO 700 

VAL8=CTHL*RADL*WINC2/HINC3 

VAL10=VAL8/2.0 



IF(I.GE.(U+1) ) GO TO 730 
VAL18=CTHM*RADH*WINC2/HINC1 
GO TO 750 

IF(I.LE.(U + ZW3)) GO TO 740 

V AL1 8=CTHL*RADH*W INC2/HI NC3 
GO TO 750 

V ALL 8=C THL*RADH*W I NC2/HI NC2 
VAL20=VAL18/2.0 



IFU.NE.UU+D+ZW3)) GO TO 760 
VAL5=(CTHL*(RADH**2-(R0+B)**2)+CTHM* ( (RO+B) 
***2-RADL**2 ) )/ (2.0+WINC1 ) 

GO TO 800 

IF(I .GT.( < U+l )+ZW3) ) GO TO 770 
V AL5 = CTHM*ALI NE/ ( 2 .0*W I NCI ) 

GO TO 800 

VAL5=CTHL*ALINE/(2 .O^WINCU 



I F ( I .GE . ( U+l ) ) GO TO 810 
VAL7 = CTHM*ALINE/( 2 .0*WINC2) 

GO TO 845 

IF (I .NE . (U + l) ) GO TO 820 

VAL7=(CTHL*(RADH**2-R0**2 ) + CTHM*( R0**2-R AOL* *2) ) 
*/ (2.0*WINC2) 

GO TO 845 

VAL7 = CTHL*ALiNE/( 2.0*WINC2) 

COUPLE CONDUCTANCE VALUES TO CONNECTION TERMS 

IF(I.EQ.l) GO TO 905 
IF(I.EQ.CO) GO TO 912 
Q = N*( I-l) 

DO 900 J=Q+1 tQ+N 



IF(J.EQ.(Q+1) ) GO 
I F ( J . LT . ( Q+ ZW + 1 ) ) 
IF (J .EQ. (Q + ZW + 1) ) 
IF (J .NE. (Q+N) ) GO 



TO 875 
GO TO 880 
GO TO 885 
TO 890 
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c 

C K ( J ♦ 1) =VAL ID 
CK( J,2)=VAL7 
CK<J,3i=VAL20 
GO TO 900 
C 
C 

880 CK(J,1)=VAL4 
CK(J,2)=VAL5 
C K ( J* 3) = VAL5 
CK< J,4)=VAL14 
GO TO 900 
C 
C 

885 C K ( J 1 1 ) =VAL6 
C K ( J» 2) = VAL5 
CK( J,3)=VAL7 
C K ( J » 4) =VAL 16 
GO TO 900 
C 
C 

890 CK(J,1)=VAL8 
CK( J,2)=VAL7 
CK(J,3)=VAL7 
C K ( J » 41 =VA L 18 
GO TO 900 
C 
C 

875 CK ( J 1 1) =VAL2 

CK ( J » 2 ) =VAL5 
C K ( J » 3j =VAL 12 
C 

900 CONTINUE 

901 GO TO 1000 
C 

C 

905 00 910 J = 1 » N 
IF(J.GT.l) GO TO 906 
CK( J, 1) =V AL5 

C K ( J » 2) =VAL12 
C K ( J » 31 =A 1 
GO TO 910 
C 
C 

906 IF (J.GE.ZW+11 GO TO 907 
CK(J»1)=VAL5 
CK(J»2J=VAL5 

C K ( J » 3 ) =V AL 14 
C K ( J ♦ 4) =A2 
GO TO 910 
C 
C 

907 IF(J ,NE.(ZW+in GO TO 90 8 
CK(J, 1)=VAL5 
CK(J»2)=VAL7 

C K ( J » 3) =VAL16 
CKI J,4J=A3 
GO TO 910 
C 
C 

908 IFIJ.EQ.NJ GO TO 909 
C K ( J ♦ 1 ) =VAL7 

CK( J,2J=VAL7 
CK( J,3)=VAL18 
C K ( J ♦ 4) =A4 
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GO TO 910 
C 
C 

909 CK ( J» 1) =VAL7 
CK ( J , 21 = VAL20 
CK1J, 3)=A5 

C 

910 CONTINUE 
C 

C 

91 1 GO TO 1000 
C 

C 

912 Q= (CO— 1 1 *N 

DO 917 J=Q+1 • Q+N 
IF (J .GT. (Q+1J ) GO TO 913 
CK! J, 11=VAL2 
C K ( J f 2 ) =V AL5 
CK ( J , 3) =A1 
GO TO 917 
C 

c 

913 IF<J.GE.lZW+3+l)) GO TO 914 
CK ( J » 1 ) =V AL4 

CK ( J» 2) =VAL5 
CK ( J, 3) = VAL5 
CK ( J,4) =A2 
GO TO 917 
C 
C 

914 IF (J.NE.l ZW+O + DJ GO TO 915 
CK ! J , 1 ) =VAL6 

CK! J,2)=VAL5 
CK ( J t 3) =VAL7 
CK ( J , 4) = A3 
GO TO 917 
C 
C 

915 IF(J .EQ. (Q+N)) GO TO 916 
C K ( J » 11 =VAL8 
CK(J,2)=VAL7 

CK ( J,3) =VAL7 
C K ( J » 4) =A4 
GO TO 917 
C 
C 

916 CK ( J, 1) = VAL 10 
CK 1 J » 2 ) =V AL 7 
CK(J»3)=A5 

917 CONTINUE 
1000 CONTINUE 
C 

C 

C END OF CONDUCTANCE CALCULATIONS 
WRITE ( I OT » 922 ) 

922 FORMAT!/,* ENTER NAME OF INPUT FILE*:) 

READ! IN, 923 )FNAME 

923 FORMAT ( A12 ) 

C 

C 

C CREATE DATA FILE FOR TVSSI 

r 

OPEN (IWR,FILE=F NAME, STATUS 3 'NEW* , FOR M= • F ORMATTE D* , 
*IOSTAT=ICK) 
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C READ IN NODE TEMPERATURES GENERATED BY TVSSI 

C 

C 

DO 924 J=1 , 25 

READ! 16,1045) ( T< I ) , 1= ( ( J- 1 ) *12+1 ) , J*1 2 J 

924 CONTINUE 

C 

WR ITE (I OT , 925 ) 

925 FORMAT!/' WRITING FILE TVSSI') 

C 

WRITc(I0T»92£» ) CO#N »C0 »N , ZW2 

926 FORMAT!/' THIS ',13,' NODE MODEL HAS ',12,* ROWS • 
*' AND* , 12 » ' COLUMNS. THERE ARE ',12,' ROWS') 

WRITEdOT , 927 ) ZW3 

92 7 FORMAT!' IN THE CONDESATE LAYER AND ',12, 

*• IN THE FIN. THE FIN CONTAINS') 

WRITE(I0T,928) ZW 
928 FORMAT! 12 » ' COLUMNS. ') 

C 

C COMPUTE HEAT FLOW ACROSS DESIGNATED SURFACES 

C 

QF L=0 .0 
QFS=0.0 
Q F U= 0 . 0 
QU=0 . 0 
Q0=0 . 0 
C 
C 

DO 2000 1 = 1, N 

949 IFII.GT.! N— ZW-1 ) ) GO TO 955 

V 1=( U + l ) *N+ i+ZW+I 
V 2 = ( U*N) +1+ZW+I 
IF1V2.EQ. (U+1)*N) GO TO 950 
QFL=OFL+( T ! VI )-T! V2 ) ) *CK ( V2 ,4) 

GO TO 955 

9 50 QFL=QFL-HT< VI )-T ( V2))*CK ( V2 ,3) 

955 IF1I.GT.ZW3) GO TO 960 
W=N*!U+I J+2+ZW 

QFS=QFS+!T(W)-T!W-1))*CK( W, 3) 

960 I F ( I • GT • ( Z W + I ) J GO TO 970 
X I=! U+ZW3 ) *N+ I 
X2=! U+ZW3+1 )*N+I 
IF1I.EQ.1) GO TO 965 
QFL= QFU + ! T ( X2 )-T t X 1 ) ) *CK ! XI , 4 ) 

WRITEI6 ,9000) QFU, T 1X2 ) , T 1 XI ) , CK< XI, 4 ) , XI , X2 
GO TO 970 

965 QFU=QFU + (T! X2 )-T( X 1 ) ) *CK ( X 1 , 3 ) 

WRITt!6,9000)QFU, T!X2) ,T ( XI ) , CK( XI, 3 ) , XI , X2 
9000 F0RMAT(1X,4F 10.5,213) 

970 I F ( I . GT . 1 ) GO TO 980 
975 QO=QO+(T< I)-TINJ*CK!I ,3) 

QU=QU+(TS-T ( I + ! N* ( CO-1 ) ) ) ) *CK ! ! I + ! N*( CO- 1 ) ) ) ,3) 

GO TO 2000 

980 IF(I.EQ.N) GO TO 975 

QO=QO+(T! I)-TIN)*CK( 1,4) 

QU=QU+! TS-T ( I + ( N* ( CO- 1 ) ) ) ) *CK ( ( I +( N* ( CO- 1 ) ) ) ,4) 
2000 CONTINUE 

QT 1=QFU+QFS 
QT0=QT1+QFL 
C WRITE STATEMENTS 
WRITEdWR, 1030) 

WRITE! IWR, 1031 )RI , RO, B 
WRITEIIWR, 1032 ) 

WRITEdWR, 1033) C,D,A 
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OOh OO 



c 

10 30 

1031 
10 32 

1033 
10 34 

1035 

1036 

10 37 

10 40 

1041 
10 42 

1043 

1044 

10 45 
1047 

10 48 



0 51 



WRITEdWR, 1034) 

WRITE II WR, 1035 )TI N,TS 
WRITEIIWR, 1036 ) 

WRITEdWR, 103 7) CTHM.CTHL , HI ,HS 
WRITEdWR, 1 047 ) QF U 
WRITEdWR, 1 048 )QF S 
WRITEdWR, 1040JQT 1 
WRITEdWR, 104DQFL 
WRITEdWR, 1042 )QT0 
WRITEdWR, 1043 )0U 
WRITE (IWR, 1044)Q0 

FORM ATI / ' PIPE INNER RADIUS PIPE OUTER RADIUS ', 
*• FIN HEIGHT ') 

F0RMAT(4X,F8.5 ,' M • ,1 IX , F 8. 5 , ' M • , 5 X , F 1 0 . 5 , • M* ) 
FORMAT!/* FIN SPACING FIN WIDTH *, 

*• CONDENSATE HEIGHT •) 

F0RMAT(1X,F8.5 ,* M • , 10X , F 8 . 5, ' M*,7X,F9.8,' M*) 

FORMAT! / * TEMPERATURE INSIDE PIPE TEMPERATURE 

*' OUTSIDE CONDENSATE LAYER •) 

F 0 RM AT ( 7X , F4 . 1 , * C ' ,34X , F 4. 1 , * C*) 

FORMAT!/* K METAL K LIQUID H INSIDE PIPE *, 

*' H IN CONDENSATE LAYER') 

FORMAT! F5.0,' W/M • ,4X , F5 . 1 , * W/ M • , 3 X , F 10 . 0, • W/M**2' , 

* 5 X, 10.0,' W/M**2') 

FORMAT!/* Q IN ON FIN UPPER AND SIDE SURFACES ', 

*F9 .3 , • WATTS •) 

FORMAT!/' Q IN ON PIPE SURFACE *,15X,F9.3,' WATTS') 
FORMAT!/* TOTAL Q IN THROUGH PIPE AND FIN SURFACE *, 
*F5 .3, ' WATTS • ) 

FORMAT!/* Q IN THROUGH CONDENSATE *,12X,F9.3,* WATTS') 
FORMAT!/* Q OUT THROUGH PIPE INNER SURFACE ',3X,F9.3, 
*• WATTS *) 

FORMAT! 12F6.2 ) 

FORMAT!/* Q IN THROUGH JPPER FIN SURFACE ', 9X , F9 . 3 , 

*• WATTS *) 

FORMAT!/* Q IN THROUGH FIN SIDE SURFACE *,9X,F9.3, 

*' WATTS •) 



CLOSE(IWR,IOSTAT= ICK) 

IF1ICK .NE. 0) WRITE! I0T,1051)EEX 
FGRMATCTROUBLE CLOSING PRINTER OUTPUT FILE* 



STOP 

END 
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